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Abstract 

The massive Schwinger model may be analysed by a perturbation expansion in the 
fermion mass. However, the results of this mass perturbation theory are sensible only 
for sufficiently small fermion mass. By performing a renormal-ordering, we arrive at a 
chiral perturbation expansion where the expansion parameter remains small even for large 
fermion mass. We use this renormal-ordered chiral perturbation theory for a computation 
of the Schwinger mass and compare our results with lattice computations. 
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1 Introduction 



The Schwinger model ]!], [|, or two-dimensional QED with one fermion flavour, is one of 
the most widely studied quantum field theoretic models. It has the attractive feature of 
being rather simple, while, nevertheless, it has many properties similar to those of more 
complicated gauge theories. Among them are the chiral anomaly 0-0, the formation 
of a chiral condensate and the non-trivial vacuum structure vacuum) [0, |[, |8|j-fT2[ . 
In addition, the fundamental fermion is confined, and only mesonic states occur in the 
physical particle spectrum [|1], [|. 

The Lagrangian density of the model is 



X 



When the fermion is massless (m = 0), the model may be solved exactly, and it turns 
out to be equivalent to the theory of one free, massive boson field ("Schwinger boson") 
with Schwinger mass /x = e/y/TT. 

For the massive (m 7^ 0) Schwinger model, the Schwinger boson turns into an interact- 
ing particle, and its mass acquires corrections. The investigation of the massive Schwinger 
model started more than 20 years ago with a series of now classical papers |]TB|-[|r7| and 
has continued since then (see |18| for a review). Quite recently, the corrections to the 
Schwinger mass, as well as the formation of higher bound states and the chiral conden- 
sate in the massive Schwinger model, have been studied by a variety of methods. Among 
these methods are mass perturbation theory ||18||— 122||, light-front quantization |23]— 128], 
lattice computations P9|- |[34| , and a generalized Hartree-Fock approach f55|— The 



multi-flavour case has been discussed, e.g., by f37|-f42|]. 

In this article we want to compute the Schwinger mass of the massive Schwinger model 
by a method that is similar to mass perturbation theory. By performing a renormal- 
ordering, we shall find a new expansion parameter (instead of the fermion mass m) 
that remains rather small even for large m. Therefore, the resulting chiral perturbation 
expansion will produce a Schwinger mass that tends to the result of mass perturbation 
theory for small fermion mass m, but remains reasonable even for large m. Finally, we 
shall compare our results to some lattice data. 

We use the bosonized, Euclidean version of the massive Schwinger model in the sequel, 
and our conventions for the massive scalar propagator are 



2tt 



K (fi\x\), D M (p) 



-1 



p 2 + jj? ' 



(2) 



where K is the McDonald function. 
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2 Renormal-ordered chiral perturbation theory 

Our starting point is the Euclidean, bosonized Lagrangian density of the massive 
Schwinger model |TB|, [T3|, |T7| , 



L = —N„ 



-$(□ - yU 2 )$ + —fim cos( 
2 2n 



47T$ 



(3) 



Here /x = ejix 1 ! 2 is the Schwinger mass of the massless (m = 0) Schwinger model, 9 is the 
vacuum angle and 7 = 0.5772 is the Euler constant. Further, denotes normal-ordering 
w.r.t. [L. 

By making use of the well-known normal-ordering relation 



fx 



(4) 



(see ||43l) 1 we may rewrite the Lagrangian density (3) for arbitrary normal-ordering mass 
p, like 



L = -Nn 



1 e 7 
-$(□ - /i 2 )$ H /2mcos( 

2 ' 27T 



47T$ 



(5) 



up to an irrelevant constant. Next, let us observe that the interaction term ~ cos\/47r$ 
contains a $ 2 piece, too. We now want to shift this $ 2 piece from the interaction term 
to the mass term in (5), and we want to choose the normal-ordering mass ft equal to the 
total mass of the resulting mass term. We find the equation 



P 2 



/j 2 + 2e 7 cos 9 jim 



with the solution 

Jx = e 7 cos 9 m + 
Our new Lagrangian density reads 



jj? + e 27 cos 2 9 m 2 . 



(6) 
(7) 



-Np [-$(□ - /i 2 )$ + Xi-^e^* + e~ w e 



+ 27T COS 9 $ — COS I 



Z71 



(9) 



where we have subtracted a constant (cos#) in the interaction term in order to avoid 
infinite vacuum energy contributions to the perturbation series. Further, we rewrote the 
sine and cosine of the field $ in terms of exponentials, because the latter ones have simpler 
Wick contractions (see (12) below). By this simple rewriting of the Lagrangian density 
we have achieved a major advantage. A perturbative expansion of the original version 
(3) of the Lagrangian is an expansion in powers of m and will, obviously, only lead to 
sensible results for sufficiently small m. On the other hand, for the new version (8) of the 
Lagrangian, the (dimensionless) expansion parameter of the perturbation expansion is 



X_ 

J 2 



e 1 m 



e 1 m 



2-kli 27r(e 7 cos 9 m + ^ li 2 + e 27 cos 2 9 m 2 ) 



(10) 
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Figure 1: The figure shows how the expansion parameter (y-axis) varies with the 
fermion mass m. We use units — 

which is rather small for all m provided that the vacuum angle 9 is not too large. E.g., 
for 9 = we find 

= (A)( m = , = 0) < 9 = 0) < (A)( m = oo,9 = 0) = J- (11) 

(see Fig. 1). 

Therefore, the perturbative expansion for the theory defined by (8) leads to finite 
results for all m and, as the expansion parameter is quite small, one may hope that these 
results are at least qualitatively correct even for large m. The interaction Lagrangian 
in (8) contains a $ 2 piece in addition to the exponentials, therefore the resulting chiral 
perturbation theory is slightly more complicated than the usual mass perturbation ex- 
pansion, because it contains Wick contractions among both polynomials and exponentials 
in $. All these may be inferred from the general formula (see [IT], pE2| ) 



n 



i=i 



■ e~\§ d,yid,y2p(yi)Df i (y 1 -y2)p(y2)+ii3j2k=i ak I d yp(y) D p(y- x k) Q2) 

where <j; = ±1 and p(x) is a source term for the <&(x) field (i.e., <$>(xi) are inserted into 
the VEV (12) by taking functional derivatives 5/5p(xi)). The subscript means A = 0, 
and (5 = \Z4tt in our case. Further, the exponentials on the l.h.s. are normal-ordered w.r.t. 
the mass of the $ field, p,. 

The Lagrangian density (8), together with the Wick contraction formula (12), will 
now serve as a starting point for the perturbative computation of the two-point function 
($(?/i)$(?/2))a and for the approximate determination of the Schwinger mass. 
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3 Two-point function and Schwinger mass 

Before actually starting the computation, let us observe that already the zeroth order 
Schwinger mass fx gives a certain prediction for large m if 9 = (numerical estimates for 
large m will always be done for 9 = 0, because for 9 ^ the expansion parameter \/jx 2 is 
larger in the large m region, making thereby numerical estimates less reliable). We find 

lim u(m, 9 = 0) = 2e 7 m ~ 3.56m. (13) 

The true value is 2m, because in the limit of vanishing (electromagnetic) interaction the 
Schwinger boson $ is just composed of two free fermions, each of mass m. Therefore, 
this lowest order estimate is nearly 80% off the exact value (which is not too bad for a 
lowest order estimate). Further, when Taylor-expanded in m, p reproduces the first order 
correction of mass perturbation theory (M . . . Schwinger mass in zeroth order) 

o 

TfL 

M := p = /i + e 7 cos 6»m + o(— ). (14) 

A* 

This feature should generalize to higher orders: The two perturbation expansions (in m 
for the Lagrangian (3); in A for the Lagrangian (8)) are just different formulations of the 
same theory; further, the parameter A (or equivalently \/p 2 ) has a well-defined Taylor 
expansion in m, starting with a linear (m 1 ) term. Therefore, the order n result of the 
perturbation series in A, when expanded in m up to n-th order, should coincide with the 
n-th order mass perturbation theory result. We conclude that there is no A 1 contribution 
to the A perturbation series for the Schwinger mass. 

The perturbation expansion for the two-point function is 

(*(i/i)*(ifa)>A = ^(HyiMv2) E -t II / ^i(^)>o, (is) 

Z n=0 U - 1=1 J 

but we prefer not to include the 1/Z normalization factor explicitly, and we will remove 
disconnected pieces by hand. Therefore we want to compute 

oo \n n „ 

(*(l/i)*(!fc)>A = (*(vi)*(lfc)>o + <*(vi)S(lfc) E -T II / dxM*i))o- (16) 

n=l n - l=l J 

First we observe that the A 1 contribution is indeed zero. This may be checked by ex- 
plicit computation, but it can be understood immediately from the following observation. 
When only one L\ is inserted into the two-point function (16), then for the exponentials 
exp ±i\^4irQ(x) in L\(x) only the $ 2 part is effective, because a higher number of boson 
lines emerging from the vertex L\ cannot be Wick-contracted to anything. But precisely 
this $ 2 piece is subtracted in Lj, (8), therefore the net contribution is zero. 
Up to second order, we obtain 

A 2 f 

(<%i)$(y 2 )>A = -D p ( yi -y 2 ) + — J <foi<fc 2 ($(yi)$(y 2 )- 
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.( e * giV^si) + e -te e -iV^Hxi) + 47r cos ^ $2( Xi ) _ 2 C o S 0). 

• ( e <« e *^*(* a ) +e -io e -iV&*(x2) +47rcos^$ 2 (a; 2 ) -2cos#), (17) 

where we now have to perform all the Wick contractions with the help of formula (12). 
Using the symmetry Dp(x) = D^—x) and the relabelling symmetry / dx\dx2f(xi,X2) = 
J dxidx2f(x2,xi) and omitting the disconnected pieces, we arrive at (c . . .connected) 

($(3/i)$(j/ 2 ))a = -D p ( yi - y 2 ) + 2rr\ 2 J dx x dx 2 - 

■ - cos 29e^ D ^-^\D-,{ yi - x 1 )D ii (y 2 - x 2 ) + D p ( yi - xjD^ - x,)) 
-e-^(xi-x2)(_ jD . (j/l _ Xl ) D -(y 2 - X2 ) + D -( yi - Xl )D p (y 2 - m)) 
+ (cos2# + - x 1 )D ii {y 2 - an) + 2^^ - y 2 )^(^i - x 2 ) 

+87rD /i (xi - x 2 )D {i (y 1 - xi)D Vb (y 2 - x 2 ) + 8tt 2 D?(xi - x 2 )D fl (y 1 - xi)D n (y 2 - an)) 

+ (cos2# + - a;i)Dp(j/ 2 - ari) 

-4tt(cos 26> + lj-D^an - x 2 )D il (y 1 - x 1 )D jl (y2 - x 2 ) 

- (cos 2^ + l)£> p (yi - Xi)Dp(j/ 2 - an)] , (18) 

where in the A 2 part the first two terms (with the exp ±4TrDp(xi — x 2 ) factors) are from 
the (<3><3>exp <I> exp $) contractions, the next term is from the ($$$ 2 exp<3>) contractions, 
and the three last terms are from the (<3><3> exp <E>), ($$$ 2 $ 2 ) and ($$<1> 2 ) contractions, 
respectively. It is important now that in expression (18) the exponentials exp rkAnD^xi — 
x 2 ) = 1 ± AttD^xi — x 2 ) + . . . that are integrated separately get their constant piece 1 
cancelled by some other contributions, and those exponentials that are convoluted with 
propagators get their first order piece iinD^xi — x 2 ) cancelled (their constant piece 
may be omitted by hand, because it is a disconnected term). Therefore, we arrive at 

m yi My 2 )) c x = -D p ( Vl - y 2 ) - 2n\ 2 - 
■ cos29(E + (fi) J dxD fl (y 1 -x)D fl (y2-x)+ J dx 1 dx 2 E^ ) {x 1 -x 2 , p)Dp,{y 1 -x 1 )Dp.{y2-x 2 )) 
+E-(fi) J dxD Pj (y 1 -x)Dp i (y 2 -x)- j dx 1 dx 2 E {1) (x 1 -x 2 , p)D p .(y 1 -x 1 )D Ji (y2-X2)) (19) 
- 8tt 2 (cos 29 + 1) J dx^D 2 ^ - x 2 )D p ( yi - x 1 )D p (y 2 - an)] , (20) 
where E±(p), E2\x,p) are defined like 

E±{x,p) = e ±4nD ^ x) - 1, E±{ji) = J d 2 xE±(x,fi) 

E ( ±\x,jl) = e ±AnD ^ x) - l+AnD-^x). (21) 
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The last term in the above expression for (§(yi)&(y 2 ))x, which we singled out by its own 
number (20), will be significant later on. 

The Fourier transform of the two-point function reads 



<*$>a(p) 



2ttA 2 



cos 20 (£+(//) + E { +\p,fi)) 



P 2 + fi 2 (p 2 + /i 2 ) 2 

+ E-{p) - E W (p,fi) - 8vr 2 (cos20 + 1)51(0) 



'.22) 



p 2 +fi 2 + 27i\ 2 (cos29(E + (fi) + E^{p, fi)) + E-(fl) - E w (p,fi) - §f (cos20 + V 

where we inverted the two-point function perturbatively and used 

1 



Dm 



r,2 



(23) 
(24) 



(-Djl(p) is just a two-boson blob). 

The Schwinger mass is given by the zero of the denominator in (23), which leads to 
the equation (after rescaling, e.g., E±\p,p) = 4^E^\f, 1), etc.) 



2ttA 2 



cos20(£+(l) +E { +\-, 1)) + £_(1) -£ w (-, 1) -2tt(cos20 + 1) (25) 
fi fi J 



H* fi 

with the second order solution (M 2 . . . second order Schwinger mass) 
M 2 2 := -p 2 \ x * = fi 2 + ^ [cos W(E+(1) + E$\i, 1)) 



+E-(1) - E W (i, 1) - 2tt(cos20 + 1) 
27rA 2 

ft 2 + —^(2ticos29A + + 2nA_ - 2tt(cos20 + 1)), 
fi 2 



where 



A 



± 



drr 



e -2K (r) _ 1 + /o ( r )( ±e T2^oW -p X + 2K (r)) 

A+ = -0.6599, A_ = 1.7277 



(26) 

(27) 
(28) 



see 



20] , pT8|| (we changed notation from A, B in |20|, [T8|] to A + , A_ here). Re-expressed in 



fi, m we find 



M\ = fi 2 + 2e 27 cos 2 9 m 2 + 2e 7 cos 9 m^J fi 2 + e 2 ^ cos 2 9 m 2 
+e 2y m 2 (A + cos 29 + A J) - e 27 m 2 (cos 20 + 1) 



/i 2 + 2e 7 cos9mJfi 2 + e 2 ^ cos 2 9 m? + e 27 m 2 (A + cos 29 + A_ ) . 



(29) 
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Observe that one piece of the second order (in A) contribution, stemming from the term 
that we singled out in (20), has precisely cancelled the o(m 2 ) piece of the lowest order 
(in A) contribution, /i 2 . This cancellation ensures that M 2 , when Taylor expanded in m 
up to second order, coincides with the second order result of mass perturbation theory 
(MPT), 

(M 2 MPT ) 2 = /i 2 + 2e 7 cos ^ /im + e 27 m 2 (A+ cos 20 + A_), (30) 



see §§ 

Numerically, when evaluated at 9 = and for large m, 

lim M 2 (m, 9 = 0) = me 7 (2 + A + + A) 1/2 ~ 3.12m, (31) 

we find that M 2 (m = oo,9 = 0) is still about 55% off the true value. Therefore, a 
higher order computation would be highly desirable. However, already the third order 
computation is very complicated. Still, we are able to compute a part of the third order 
contribution with considerably less effort. To achieve this aim, we use the requirement 
that the perturbation series in A, when Taylor expanded in m, has to coincide with the 
mass perturbation series up to the given order. The lowest order contribution (in A) to 
M 2 , p 2 , contains a third order (in m) piece, 



3 

£ 2 L 3 = e 3 -y C os 3 e— , (32) 
fi 



and such a term cannot be produced by ordinary mass perturbation theory. Therefore, 
the third order (in A) contribution must contain a term that, when Taylor expanded in 
m, precisely cancels the term (32). The contribution that does the job reads 

8 1 M 2 = -(2tt) 3 cos 3 9 — A = -e 37 cos 3 9 — . (33) 

/x 4 ft 

A further contribution of order A 3 may be identified by the same argument via a closer 
inspection of ordinary mass perturbation theory in order m 3 . For this, observe that one 
contribution to our second order result (26) stemmed from the evaluation of E±(J=, 1) at 
(p/fi) 2 = —1. For second order mass perturbation theory an identical term occurs, but 
this time evaluated at (p/fi) 2 = — 1. In addition, in mass perturbation theory there is a 
first order contribution <5(M{ vIPT ) 2 = 2e 7 cos 9fj,m, therefore there must be a third order 
contribution (in mass perturbation theory!) 

5(Mr T ) 2 ^y(cos2^«(p, 1) - EW(p, l))^ (34) 

or, including all coefficients, 

-(2e 7 cos^m)^^ (cos 29 E^(p, 1) - E^\p, 1))|„ 2= _ 1 
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e 3 " 7 cos 9 m 3 d 

7T/i dp 2 



J dxe ipx {cos26E { +\x, 1) - £ W (x, l))|p2=_i 



= 2e 37 cosfl — (cos29B + - B_) (35) 

B±= J drr 3 I (r)E ( 2\r,l) (36) 

5+ = 1.449, = 1.912. (37) 

This third order contribution must be matched by a corresponding third order term 
in A perturbation expansion (observe that this term in A perturbation theory cannot be 
produced in the same way, because there is no first order correction SMi in A perturbation 
theory; most likely, the term is produced by the ($$$ 2 exp $ exp $) contributions). This 
corresponding term reads 

\ 3 ™3 

8 2 Ml = 2 (2tt) 3 cos 0—(cos26B + - B_) = 2 cos 9e 3 ~<— (cos20£ + - £_). (38) 

Therefore, our new, partially third order corrected expression for the Schwinger mass is 

M 2 2 i = M\ + biM\ + 6 2 Ml (39) 

Numerically, for 9 = 0, M 2 i has the large m limit 

jirn^ M 2 1 (m, 9 = 0) = 2.58m, (40) 

which is nearly 30% above the exact value 2m. 

Of course, there are still some third order contributions missing, but as they are due 
to higher-dimensional convolution integrals of exponentially decaying functions, we may 
hope that their contribution is not too large. 

We display the results of our renormal-ordered chiral perturbation theory (for 9 = 0) 
in Fig. 2, and compare them with the lattice data of |3D| in Fig. 3. Fig. 3 shows 
that our result M 2 i agrees with the lattice data very well up to m ~ 1; for larger m the 
deviation increases, reaching the above-mentioned 30% for m — > oo. 

For an overview on additional data for the Schwinger mass (e.g., from light-front 
computations) we refer to |25j (especially their Fig. 2 and Table 16), and to [EH. 



4 Summary 

By a simple renormal-ordering of the bosonic Lagrangian of the massive Schwinger model 
we were able to find a new perturbation expansion parameter, (see (10), Fig. 1), that 
remains rather small even for large fermion mass m. The resulting chiral perturbation 



theory is similar to mass perturbation theory pOl 18 1, although it involves somewhat 



more tedious computations. We used this renormal-ordered chiral perturbation theory to 
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compute the Schwinger mass M, including all second order (in ^) and some third order 
contributions. The resulting expression for the Schwinger mass M reduces to the result 
from mass perturbation theory (||2U||) when Taylor-expanded in the fermion mass m. For 
large m our result describes the known data for the Schwinger mass (e.g. the lattice data 
from qualitatively correct, although there remains some deviation (e.g. our result 
M 2 i differs by nearly 30% from the known value 2m at the far end m — > oo). 

Actually, we may conclude from Fig. 3 that high order (in ^) perturbative contribu- 
tions will be relevant for m — > oo. The reason for this is that the expansion parameter 
is nearly constant for ra 6 [1, oo] , see Fig. 1. Therefore, low order contributions to M will 
shift the curves for M 2 ± (or M 2 ) in Fig. 3 by a nearly constant amount for m E [1, oo]. 



1 than at m 



oo, 



However, the curves in Fig. 3 are closer to the lattice data at m 

and they can only approach the lattice data everywhere when high orders in -fp give a 
non-negligible contribution. 

Finally, it should be mentioned that it could be interesting to apply this renormal- 
ordered chiral perturbation theory to the computation of other physical quantities of the 
massive Schwinger model (e.g., higher bound states, chiral condensate). In addition, it 
might be possible to use a version of this renormal-ordered chiral perturbation theory 
for the massive multi-flavour Schwinger model, where ordinary mass perturbation theory 



fails due to infrared divergencies 41, 42 
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Figure 2: The figure shows successive perturbative results for the Schwinger mass, fi is 
the upper curve, M 2 is the curve in the middle and M 2 i is the lowest curve. We use the 

rescaled Schwinger mass M' = M/(l + m 2 ) 1 / 2 , like in [25, |26|, and choose units \i = 1. 
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Figure 3: The figure compares our results M 2 , M 2 i to the lattice data of [^] (small dots, 

We again use the 



for larger m) and (fat dots, for small m, with rather large errors) 
rescaled Schwinger mass M' = M/(l + m 2 ) 1//2 and units fj, — 1. 



